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Abstract 

Invasion fronts in ecology are well studied but very few mathematical results concern the case with variable 

motility (possibly due to mutations). Based on an apparently simple reaction-diffusion equation, we explain the 
observed phenomena of front acceleration (when the motility is unbounded) as well as other quantitative results, 
such as the selection of the most motile individuals (when the motility is bounded). The key argument for the 
construction and analysis of traveling fronts is the derivation of the dispersion relation linking the speed of the 
wave and the spatial decay. When the motility is unbounded we show that the position of the front scales as t^^^. 
When the mutation rate is low wc show that the canonical equation for the dynamics of the fittest trait should 
be stated as a PDE in our context. It turns out to be a type of Burgers equation with source term. 

Resume 

Fronts d'invasion avec motilite variable : repartition des phenotypes et acceleration de I'onde. Les 

fronts d'invasion en ccologic ont etc largomont otudics. Ccpcndant pcu de rcsultats mathcmatiqucs existent pour 
le cas d'un coefficient dc motilite variable (a cause des mutations). A partir d'un modele minimal de reaction- 
diffusion, nous expliquons le phenomene observe d'acccleration du front (lorsque la motilite n'est pas bornee), et 
nous demontrons I'existence d'ondes progressives ainsi que la selection des individus les plus motiles (lorsque la 
motilite est bornee). Le point cle pour la construction des fronts est la relation de dispersion qui relie la vitesse de 
I'ondc avec la dccroissancc on ospaco. Lorsque la motilite n'est pas bornee nous montrons que la position du front 
suit une loi d'echelle en t^^^. Lorsque le taux de mutation est faible, nous montrons que, dans notre contexte, 
I'equation canonique pour la dynamique du meilleur trait est une EDP. C'est une equation de type Burgers avec 
terme source. 
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Version frangaise abregee 



Dans cette note nous etudions un modele simple decrivant des fronts invasifs en ecologie, pour lesquels 
la mobilite des individus est sujette a variations. Le modele, issu de [1], est le suivant, 

{dt'n{t, X, 6) — 9dl^n(t, x, 9) + rn(i, x, 9) (1 — p{t, x)) + adggn{t, x,9) , x € (— oo, +oo), ^ > , 
p{t,x) = J n{t,x,6)de , n{t,-oo,9) =N{9) , n{t, +oo, 9) = . 

Les conditions aux limites sont completees ci-dessous. Nous etudions la dynamique de cette equation sous 
differents regimes. Dans un premier temps nous etudions le probleme de propagation de front a vitesse 
constantc. Ccla circonscrit I'cspacc des traits 9 G (0, 0), qui doit ctrc borne O < +oo. Nous dcmontrons 
que la situation est similaire au cas de I'equation de Fisher-KPP. II existe une vitesse minimale c* de 
propagation des ondes de reaction-diffusion (Theoreme 1). La relation de dispersion, analogue a celle de 
Fisher-KPP (un trinomc du second dcgrc cn roccurrcncc) est donncc via la resolution d'un problcmc 
spectral. On lit sur la distribution des phenotypes (le vecteur propre) que les fortes motilites sont favo- 
risees, conclusion opposee au cas des domaines bornes en espace [2]. Ce meme probleme spectral intervient 
lorsqu'on ctudic la propagation du front cn regime asymptotique hyperbolique (t, x) — > (t/e, x/s) dans la 
limite WKB de I'optique geometrique, en suivant I'ansatz n^{t,x,9) = exTp{u^{t,x)/e)N'^{t,x,9). 

Dans un second temps, nous considerons un espace des traits non borne, 9 = +oo. Dans ce cas 
le front accclcre sans cesse et nous montrons hcuristiqucmcnt que la loi de propagation du front est 
naturellement {x) ~ (a^/^r'^/^) t^^'^. Nous exhibons une solution particuliere qui confirme cette heuristique 
(Proposition 2). 

Dans un troisicmc temps, nous etudions le regime de mutations rarcs, et nous ccrivons une equation 
canonique pour revolution du trait selectionne localement a I'avant du front en regime asymptotique. La 
derivation formelle de cette equation conduit a une equation de transport de type Burgers, avec terme 
source (Proposition 3) . La partic transport est diic a la progression du front qui " transportc" les individus 
et done le trait selectionne localement. Le terme source est du a la pression de selection qui tend a faire 
augmenter la valeur du trait selectionne localement. 



Introduction 

Recently, several works have addressed the issue of front invasion in ecology, where the motility of 
individuals is subject to variability [3,4]. It has been postulated that selection of more motile individuals 
can occur, even if they have no advantage regarding their reproductive rate (spatial sorting) [5,6,7,8]. This 
phenomenon has been described in the invasion of cane toads in northern Australia [9] . It has been shown 
that the speed of the front increases, coincidentally with significant changes in toads morphology. Up 
to now, only numerical simulations have been proposed to address this issue. Here we analyse a simple 
model of Bcnichou et al [1] which contains the basic features of this process: spatial mobility, logistic 
reproduction, and variable motility. It is given by equation (1) where the last term on the right hand 
side accounts for modifications of the dispersal rate 6 of individuals due to mutations. We consider that 
mutations are random and they act as a diffusion process in the phenotype space. When needed, we 
impose Neumann boundary conditions in the variable 9 and far-field conditions in the variable x. 

Email addresses: emeric.bouiniaens-lyon.fr (Emeric Bouin), vincent.calvezaens-lyon.fr (Vincent Calvez), 
nicolas.meimierSparisdescartes.fr (Nicolas Meunier), mirrahimiScmap. polytechnique.fr (Sepideh Mirrahimi) , 
benoit.perthame8upmc.fr (Benoit Perthame), g . raoulSdamtp . cam . ac . uk (Gael Raoul), voiturie@lptmc.jussieu.fr 
(Raphael Voituriez). 
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1. Phenotype selection and spatial sorting in the traveling wave 



We first consider bounded dispersal rates, say 9 G (0, 9). Following [1], we seek a traveling wave solution 
to the equation (1) connecting to the uniform stationary state N{9) = 6~^. For x large, we make the 
following ansatz n{t,x,9) = cxp(A(a; — ct))Q{9), where c > is the speed of the wave, A < is the 
spatial decay and Q{9) denotes the phenotypic distribution of the individuals at the edge of the front. 
The dispersion relation is equivalent to the following spectral problem: Given a spatial decay rate A < 0, 
find c(A) and an eigenvector Q{9,X) such that 



(Ac(A) + ^A^ + r) Q{e, A) + adl0Q{9, A) = , 
deQ{0,X) = deQ{e,X) = 0, y9Q{9,X)>0, 



J QiO,X) 



d9 = l. 



(2) 



The invasion front speed c(A) is such that is the principal eigenvalue of this spectral problem. Like the 
Fisher-KPP equation, there is a minimal speed c* > associated with a critical spatial decay A* < 0. 

Theorem 1 (Front propagation and spatial sorting) For all c > c* , there exists formally a trav- 
eling front solution n(t,x,9) — N(x — ct,9). It satisfies N(z,9) ^ cxp{Xz)Q{9, X) as z ^ +oo. The 
phenotypic distribution at the edge of the front Q{9, A) is unbalanced towards more motile individuals. 
More precisely, we have 

(^)edge(A) := J 9Q{9, X)d9>^, (0),,^, (A) -^^ & . (3) 

This problem is closely related to front propagation in kinetic equations [10]. There is a natural extension 
of this result for other mutation operators as integral operators a J G{9, 9')n{t, x, 9') d9' . In this case, the 
solution of the spectral problem is deduced from the Krein-Rutman theorem. 

Interestingly, we can measure the asymmetry of the phenotypic distribution Q{9,X). The relevant 
quantity here is the mean diffusion coefficient at the edge of the front (0)edge(A). In order to show (3), we 
integrate the spectral problem (2) over (0, 0). We get Ac + A^(^)edge + r = 0. Dividing by Q{9, A), and 
integrating again the spectral problem we get after integration by parts, 



(Ac + r)e + ^A2 



+ a 



dgQ{9,X) 



Hence, the mean diffusion coefficient is given by 

9 



(^)edge(A) 



+ 



a 



9A2 



/ 



Q{9,X) 



deQ{9,X) 



d9 = 0. 



Q{9,X) 



dO > 



9 



(4) 



Individuals at the edge of the front are more motile than at the back of the front. There, the population is 
homogeneous since 7V(— oo, 9) = and the averaged diffusion coefficient is (6')back = 9/2. The estimate 
(4) measures how far the phenotypic distribution differs from the uniform distribution. Finally we find 
easily that as a vanishes, the distribution Q{6, A) concentrates as a Dirac mass on ^ = 9 hence the 
second statement in (3): the most motile individuals are selected. The conclusion is exactly the opposite 
on bounded spatial domains [2]. 

We can derive additional information for the minimal speed c*. Differentiating (2) we obtain 



(Ac'(A) + c(A) + 26IA) Q{9, A) + (Ac(A) + 9X^ + r] 



—Q{9,X)+adlg—Q{9,X) 



0. 



Definition of c* ensures that c'(A*) = 0. We use the notation (/) := J f{9)Q* (9) dO. Since the operator 
in (2) is self-adjoint, we obtain after multiplication by Q{9, A) and integration, the relation 
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c*{Q*)+2\*{9Q*)=Q, c* = -2X*^^, (5) 
Recalling that A*c* + (A*)^(^) + r = 0, we can eliminate A* and wc get the following expression for c*, 

(c.,^=4.w(l-(l-<|M)y">4.W^ 

In other words, the usual formula for the KPP wave speed underestimates the actual minimal speed. 



2. Spatial sorting and the invasion front 

Next, we focus on the invasion front. It is natural to perform the hyperbolic rescaling {t, x) {t/e, x/e) 
in order to catch the motion of the front [11,12]. The new equation writes after rescaling 

edtrf{t, X, 6) = e^Odl^n^t, x, 6) + rrf{t, x, 6) (1 - x)) + adl^rf{t, x, 6) . 

Wc perform the partial WKB ansatz in the x variable: 0) = CKY>{v^{t,x)/e)N'^{t,x,0), with the 

renormalization J N'^(t, .t, 9) dd = I. As e ^ 0, the first order expansion leads to solve 

dtu^it, x)N" = 6*19^.1*" (f, x)\^N"{t, X, 9) + r {l - p°{t, x)) N^{t, x, 9) + ad^gN^{t, x, 9) . 

The edge of the front is delimited by the area {vP{t,x) < 0}. On this set wc have p^{t,x) = by 
construction. Therefore we shall solve again the spectral problem (2) for iV° > 0. Consequently the 
motion of the front is driven by the eikonal equation built on the effective speed c(A), 

max (u", dtu^ + (9^M° • c{d^vP)) = . 

The rigorous derivation of this Hamilton- Jacobi equation requires some work. We need basically refined 
a priori estimates on (m^)^. We formally show the main argument leading to establish the viscosity limit 

of if in the set {vP{t,x) < 0}. Wc leave the complete proof for future perspectives. Let be a 
test function such that uP — v° has a strict maximum at {t'^,x°). The function — v'^ has a maximum 
at (i^,a;^), with close to Plugging into the equation satisfied by {u^,N^), namely 

[dtu%t,x) - 9\d^u^it,x)\'^ -e9dl^u%t,x)-r + Oie)] N%t,x,9) = ad'^gN%t,x,9) , 

we obtain at {t^,x^): 

[dtv°{t^, x^) - 9\d^v°{t^, x^)\^ -r + 0{e)] N%t^, x', 6) < ad%N^{f, a;^ 6) . 

Therefore, is a non-negative, non-trivial subsolution of the spectral problem (2). From the character- 
ization of the principal eigenvalue, we have 

dtv''{f,x') + d^v\f,x') ■ c{d^v^{f,x')) + 0{e) < 0. 

Passing to the limit e — > 0, we obtain that satisfies dtv^ + dxV^ ■ c{dxV^) < at {t^,x^). Therefore 
vP is a viscosity sub-solution of the eikonal equation dtvP + dxvP ■ c{dxvP) = in the interior of the set 
{vP{t,x) < 0}. The same argument shows that it is also a supersolution and thus a viscosity solution. 



3. Front acceleration 

In the case where the set of dispersal rates is unbounded, say 9 £ (0, +oo), then we cannot solve the 
spectral problem (2). There is no intrinsic speed of propagation, and the front is accelerating as time goes 
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Phenotypic distribution at the 




-1/(9 a) (2 sqrt(r «) t-f!) (2 (!+2 sqrt(r a) tr = 



(a) 




(b) 



Figure 1. (a) Phenotype selection. At the edge of the front the phenotypic distribution is not uniform: more motile individuals 
are selected. We solved numerically the relation dispersion for (r, a, 0) = (1, 1, 1). (b) Front acceleration. The set = 0} 
is plotted in the phase space (x, 0), for successive times. The far-right point X(,dgc(i) determines the location of the front. 



on. Heuristically, we expect the averaged diffusion coefficient to grow linearly with time (9) ^ (Voir) t 
(as for the Fisher-KPP equation set in the phenotype space). Hence the invasion front should scale as 
(x) ^ (a^/^r'^/*) t^/^ since the speed is given by c ~ {9)r. Therefore we perforin the asymptotic scaling 
(t, X, 9) — ?> (t/e, xje^l"^ ., 9/e) , in order to catch the motion of the front asymptotically. The equation writes 
after rescaling 

edtn'' = e^9dl^n^ + rn'' (1 - x)) + e^adlgu" . (6) 

We perform the WKB ansatz in both variables (x, 9): n^{t, x, 9) — exp x, 9)/£). We derive formally 
the Hamilton- Jacobi equation for in the limit e — > 0: 



dtu° - 9\d,u°\^ - a|aeu"|2 = r (l - x)) 



(7) 



where p^{t^ x) is the Lagrange multiplier associated with the constraint maxg u^{t, x, 9) < 0. In particular, 
when the constraint is inactive, i.e. maxg u'^{t, x, 9) < 0, we have p'^{t, x) = 0. The behavior of the phase 
function vP{t, x, 9) is as follows: the nuUset {u^{t, x, 9) — 0} propagates in space and phenotype under the 
action of space mobility, and mutations, combined with growth of individuals. We are able to compute 
explicitely the evolution of this set for a particular initial data, using Lagrangian Calculus. 



Proposition 2 (Front acceleration) For the particular initial data, u'^{0,x, 9) 



if (x, 9) ^(0,0) 
if {x, 9) ^(0,0) 



the location of the nullset {u'^{t,x,9) — 0} is given by the following implicit formula (see Fig. 1) 



9a 



((2VfS) t - 9) {29 + (2Vf^) ty 



Sketch of proof . The Hamiltonian is given by H{{x, 9), (pxjPe)) = ()\Px\'^ -\-o:\pe\'^ +r, and the corresponding 
Lagrangian writes L{{x,9), (vx,vg)) — v'^/{A9) + Ug/(4a) — r. The system of characteristics is given by 
X{t) — {29{t)p.j;{t),2apg{t)), and P{t) — — (0, |p^(t)p). Using Lagrangian formulation, we deduce after 
some calculation that, with Z the solution to equation + {129/a)Z + 24x/a — 0, x, 9) is given by 



u\t,x,9) = -±-^(9 
This enables to compute the nullset {u^{t,x,9) — 0}. □ 
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The far-right point of the curve is attained for ^edge = (\/^) t- This determines the location of the 
front. Hence the position of the front in space is exactly, 

4. Adaptive dynamics at the edge of the front 

Interestingly enough, the same equation (7) can be derived in the context of adaptive dynamics, a 
theory which studies mutation-selection processes. It is generally assumed that the mutation process is so 
slow that mutants can replace the resident species before new mutants arise, if they are better adapted to 
their environment. This yields a canonical equation which gives the dynamical evolution of the selected 
trait in the population [13,14]. Recently, PDE-based methods have been successfully used to derive such 
a canonical equation for continuous mutation-selection processes [15,16,17]. Here we extend this theory 
to the case of front invasion coupled with a basic mutation process. 

The only difference with (1) is that mutations are assumed to be rare {a — > e^a): 

dtn{t, X, 9) = edl^n{t, X, 9) + rn(i, x, 9) (1 — p{t, x)) + e'^adggn{t^ x, 9) . 

It is natural to perform a long time rescaling t — >■ t/e at the scale of evolutionary changes. Then it is 
useful to rescale space accordingly x — > a;/£ in order to catch the motion of the front (otherwise it would 
travel at order 0(l/e)). With these changes of scales, wc end up again with equation (6), resp. (7) in 
the WKB limit. We restrict to the edge of the front, namely supg vP{t, x, 0) < 0, x) = 0. We seek a 
canonical equation for the locally selected trait 6{t, x) such that u^{t, x, 6{t, x)) = supg vP{t, x, 6). 

Proposition 3 (Formal derivation of the canonical equation) The locally selected trait 6{t, x) for- 
mally satisfies a Burgers type equation with a source term, 

dtOit, x) - 2 {e{t, x)d,u°) dj{t, x) = . (8) 

The speed of the transport equation is —29{t,x)dxvP ■ It coincides with the local minimal speed of the 
traveling front, see e.g. (5). The positive source term accounts for the evolutionary drift which pushes the 
population towards higher motility (numerical simulations not shown). This equation may yield shock 
wave singularities, as for the classical Burgers equation, because more motile populations, when located 
behind less motile populations, will invade them. 

Proof. Wc start from the first order condition doU^{t, x, 6{t, x)) = 0. We differentiate this relation with 
respect to t and x, respectively, 

df^it, X, 9{t, x)) + {d'igu") dt9it, x) = , dl^it, X, 9it, x)) + (8^^) dj{t, x) = . 

On the other hand, we differentiate equation (7) with respect to 9, 

dlu° - - 29dxu°d^y - 2adeu°d^gu° = 0. 

Evaluating the latter at 9 = 9{t, x) yields 

d^tu° - 29{t,x)d^u°d^^u° = \d^u°\^ . 

Combining these calculations, we conclude that 9 satisfies equation (8). □ 
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